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About half a million years after the Big Bang, the ever-fading cosmic blackbody radiation 
cooled below 3000 K and shifted first into the infrared and then into the radio, and the 
smooth baryonic plasma that filled the Universe became neutral. The Universe then entered 
a "dark age" which persisted until the first cosmic structures collapsed into gravitationally- 
bound systems, and evolved into stars, galaxies, and black holes that lit up the Universe 
again. Some time between redshift of 7 and 15, stars within protogalaxies created the 
first heavy elements; these systems, together perhaps with an early population of quasars, 
generated the ultraviolet radiation that reheated and reionized the cosmos. The history 
of the Universe during and soon after these crucial formative stages is recorded in the all- 
pervading intergalactic medium (IGM), which is believed to contain most of the ordinary 
baryonic material left over from the Big Bang. Throughout the epoch of structure formation, 
the IGM becomes clumpy and acquires peculiar motions under the influence of gravity, and 
acts as a source for the gas that gets accreted, cools, and forms stars within galaxies, and as 
a sink for the metal enriched material, energy, and radiation which they eject. Observations 
of absorption lines in quasar spectra at redshifts up to 5 provide invaluable insight into the 
chemical composition of the IGM and primordial density fluctuation spectrum of some of the 
earliest formed cosmological structures, as well as of the ultraviolet background radiation 
that ionizes them. 
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COSMOLOGICAL REIONIZATION 



At epochs corresponding to z ~ 1000, the IGM is expected to recombine and remain neutral 
until sources of radiation develop that are capable of reionizing it. The detection of trans- 
mitted flux shortward of the Lya wavelength in the spectra of sources at 2; ~ 5 implies that 
the hydrogen component of this IGM was ionized at even higher redshifts. There is some 
evidence that the double reionization of helium may have occurred later, but this is still 
controversial. It appears then that substantial sources of ultraviolet photons were already 
present when the Universe was less than 7% of its current age, perhaps quasars and/or young 
star-forming galaxies: an episode of pre-galactic star formation may provide a possible ex- 
planation for the widespread existence of heavy elements (like carbon, oxygen, and silicon) 
in the IGM, while the integrated radiation emitted from quasars is likely responsible for the 
reionization of intergalactic helium. Establishing the epoch of reionization and reheating is 
crucial for determining its impact on several key cosmological issues, from the role reion- 
ization plays in allowing protogalactic objects to cool and make stars, to determining the 
small-scale structure in the temperature fluctuations of the cosmic background radiation. 
Conversely, probing the reionization epoch may provide a means for constraining compet- 
ing models for the formation of cosmic structures, and of detecting the onset of the first 
generation of stars, galaxies, and black holes in the Universe. 

INTERGALACTIC HYDROGEN DENSITY 

The proper mean density of hydrogen nuclei at redshift z may be expressed in standard cosmo- 
logical terms as: 



where Y is the primordial He abundance by mass, pcrit = ^Hl/iSnG) is the critical density, 
^fo = Pb/ Pcrit is the current baryonic density parameter, and Hq = 100 /i kms^^Mpc"^ is the 
present-day Hubble constant. Standard nucleosynthesis models together with recent observations 
of Deuterium yield Y = 0.247 ± 0.02, and Q^/i^ = 0.0193 ± 0.0014. Thus, 
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As some of the baryons had already collapsed into galaxies at z = 2 — 5, the value of Q.},h? 
— 0.019 should strictly be considered as an upper limit to the intergalactic density parameter. 
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Because of the overwhelming abundance of hydrogen, the ionization of this element is of 
great importance for determining the physical state of the IGM. Popular cosmological models 
predict that most of the intergalactic hydrogen was reionized by the first generation of stars 
or quasars at z = 7 — 15. The case that has received the most theoretical studies is one where 
hydrogen is ionized by the absorption of photons, + 7 — > p + e (as opposite to collisional 
ionization H + e — >p + e + e) shortward of 912 A; that is, with energies exceeding 13.6 eV, the 
energy of the Lyman edge. The process of reionization began as individual sources started 
to generate expanding H n regions in the surrounding IGM; throughout an H 11 region, H is 
ionized and He is either singly or doubly ionized. As more and more sources of ultraviolet 
radiation switched on, the ionized volume grew in size. The reionization ended when the 
cosmological H 11 regions overlapped and filled the intergalactic space. 



PHOTOIONIZATION EQUILIBRIUM 

At every point in a optically thin, pure hydrogen medium of neutral density nni, the photoion- 
ization rate per unit volume is 

nm / 7 dv, (4) 

Jul h,pU 

where Jjy is the mean intensity of the ionizing radiation (in energy units per unit area, time, 
solid angle, and frequency interval), and hp is the Planck constant. The photoionization cross- 
section for hydrogen in the ground state by photons with energy hpu (above the threshold 
hpi^L = 13.6 eV) can be usefully approximated by 

<Tn{y) = UL {yjvj:)-^, UL = 6.3 X 10-^« cm^ (5) 

At equilibrium, this is balanced by the rate of radiative recombinations p + e — > + 7 per unit 
volume, 

UeUpaAiT), (6) 

where Ug and Up are the number densities of electrons and protons, and ua — Y.io'nVe) is the 
radiative recombination coefficient, i.e. the product of the electron capture cross-section cr„ and 
the electron velocity averaged over a thermal distribution and summed over all atomic levels 
n. At the commonly encountered gas temperature of 10"^ K, = 4.2 x 10"^^ cm'^ s^^. 

Consider, as an illustrative example, a point in an intergalactic H 11 region at (say) z = 6, 
with density uh = (1-6 x 10~^ cm"^)(l + z)"^ = 5.5 x 10^^ cm^'^. The H 11 region surrounds 
a putative quasar with specific luminosity = 10"^° (z/^/z/)^ crgss^^Hz^\ and the point in 
question is at a distance of r = 3Mpc from the quasar. To a first approximation, the mean 
intensity is simply the radiation emitted by the quasar reduced by geometrical dilution, 

4-^" = 4^- m 



-4- 



We then have for the photoionization timescale: 




'CXD 



hpu 



du 



= 5 X 10'^ s, 



(8) 



and for the recombination timescale: 



1 




(9) 
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As in photoionization equilibrium nni/^ion = np/trec, these values imply nui/rip ~ 10 ^, that is, 
hydrogen is very nearly completely ionized. 

A source radiating ultraviolet photons at a finite rate cannot ionize an infinite region of space, 
and therefore there must be an outer edge to the ionized volume (this is true unless, of course, 
there is a population of UV emitters and all individual H ii regions have already overlapped). 
One fundamental characteristic of the problem is the very small value of the mean free path 
for an ionizing photon if the hydrogen is neutral, {crLnii)~^ — 0.9 kpc at threshold, much 
smaller than the radius of the ionized region. If the source spectrum is steep enough that 
little energy is carried out by more penetrating, soft X-ray photons, we have one nearly 
completely ionized H ii region, sepated from the outer neutral IGM by a thin transition 
layer or 'ionization-front'. The inhomogeneity of the IGM is of primary importance for 
understanding the ionization history of the Universe, as denser gas recombines faster and is 
therefore ionized at later times than the tenuous gas in underdense regions. An approximate 
way to study the effect of inhomogeneity is to write the rate of recombinations as 



(assuming T is constant in space), where the brackets are the space average of the product 
of the local proton and electron number densities, and the factor C > 1 takes into account 
the degree of dumpiness of the IGM. If ionized gas with electron density Ue density filled 
uniformly a fraction 1/C of the available volume, the rest being empty space, the mean 
square density would be (nl) — n^/C — {ueYC. 

The IGM is completely reionized when one ionizing photon has been emitted for each 
H atom by the radiation sources, and when the rate of emission of UV photons per unit 
(comoving) volume balances the radiative recombination rate, so that hydrogen atoms are 
photoionized faster than they can recombine. The complete reionization of the Universe 
manifests itself in the absence of an absorption trough in the spectra of galaxies and quasars 
at high redshifts. If the IGM along the line-of-sight to a distant source were neutral, the 
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resonant scattering at the wavelength of the Lya (2p — > ls]hpVa — 10.2 eV) transition 
of atomic hydrogen would remove all photons blueward of Lya off the line-of-sight. For 
any reasonable density of the IGM, the scattering optical depth is so large that detectable 
absorption will be produced by relatively small column (or surface) densities of intergalactic 
neutral hydrogen. 

GUNN-PETERSON EFFECT 

Consider radiation emitted at some frequency z/g that lies blueward of Lya by a source at redshift 
Ze, and observed at Earth at frequency z/q = ^e(l + Ze)^^- At a redshift z such that {I + z) = 
(1 + Ze)i'a/^e: the emitted photons pass through the local Lya resonance as they propagates 
towards us through a smoothly distributed sea of neutral hydrogen atoms, and are scattered off 
the line-of-sight with a cross-section (neglecting stimulated emission) of 



where / = 0.4162 is the upward oscillator strength for the transition, is the line profile function 
[with normalization / (j){u)du = 1], c is the speed of light, and e and mg are the electron charge 
and mass, respectively. The total optical depth for resonant scattering at the observed frequency 
is given by the line integral of this cross-section times the neutral hydrogen proper density nui{z), 



where di/dz = cHq^{1 + z)"^[(]m(1 + z)^ + Qk{1 + z)"^ + Qa]'^^^ is the proper line element 
in a Friedmann-Robertson-Walker metric, and ^Im, ^a, and fix = 1 — — 3re the matter, 
vacuum, and curvature contribution to the present density parameter. As the scattering cross- 
section is sharply peaked around Va, we can write 



a[iyo{l + z)] 



rUeC 



f(t>[Vo{l + z)l 



(11) 




(12) 




(13) 



In an Einstein-de Sitter (Qm = 1, = 0) Universe, this becomes 



tgp(^) 



rrieH^Vo. (1 + zfl^ 
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6.6 X lO^/i 
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nm 3/2_ (14) 

Vo.oigy ' ^ ' 



The same expression for the opacity is also valid in the case of optically thin (to Lya scattering) 
discrete clouds as long as nm is replaced with the average neutral density of individual clouds 
times their volume filling factor. 
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In an expanding Universe homogeneously filled with neutral hydrogen, the above equations 
apply to all parts of the source spectrum to the blue of Lya . An absorption trough should 
then be detected in the level of the rest-frame UV continuum of the quasar; this is the so- 
called "Gunn-Peterson effect" . Between the discrete absorption lines of the Lya forest clouds, 
quasar spectra do not show a pronounced Gunn-Peterson absorption trough. The current 
upper limit at ~ 5 is tqp < 0.1 in the region of minimum opacity, implying from equation 
(0) a neutral fraction of nHi/^n < 10~^ h. Even if 99% of all the cosmic baryons fragment 
at these epochs into structures that can be identified with quasar absorption systems, with 
only 1% remaining in a smoothly distributed component, the implication is a diffuse IGM 
which is ionized to better than 1 part in 10*^. 

In modern interpretations of the IGM, it is difficult to use the Gunn-Peterson effect 
to quantify the amount of ionizing radiation that is necessary to keep the neutral hydrogen 
absorption below the detection limits. This is because, in hierarchical clustering scenarios 
for the formation of cosmic structures (the Cold Dark Matter model being the most studied 
example), the accumulation of matter in overdense regions under the influence of gravity 
reduces the optical depth for Lya scattering considerably below the average in most of the 
volume of the Universe, and regions of minimum opacity occur in the most underdense areas 
(expanding 'cosmic minivoids'). 

A CLUMPY IGM 

Owing to the non-linear collapse of cosmic structures, the IGM is well known to be highly 
inhomogeneous. The discrete gaseous systems detected in absorption in the spectra of high- 
redshift quasars blueward of the Lya emission line are assigned different names based on the 
appearance of their absorption features (see Figure 1). 

The term "Lya forest" is used to denote the plethora of narrow absorption lines whose 
measured equivalent widths imply H i column densities ranging from 10^^ cm^^ down to 
lO^^cm"^. These systems, observed to evolve rapidly with redshift between 2 < z < A, 
have traditionally been interpreted as intergalactic gas clouds associated with the era of 
baryonic infall and galaxy formation, photoionized (to less than a neutral atom in 10^) and 
photoheated (to temperatures close to 20,000 K) by a ultraviolet background close to the one 
inferred from the integrated contribution from quasars. Recent spectra at high resolution and 
high signal-to-noise obtained with the Keck telescope have shown that most Lya forest clouds 
at 2; ~ 3 down to the detection limit of the data have undergone some chemical enrichment, 
as evidenced by weak, but measurable C iv lines. The typical inferred metallicities range 
from 0.3% to 1% of solar values, subject to uncertainties of photoionization models. Clearly, 
these metals were produced in stars that formed in a denser environmemt; the metal-enriched 
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gas was then expelled from the regions of star formation into the IGM. 

An intervening absorber at redshift z having a neutral hydrogen column density ex- 
ceeding 2 X 10^''cm~^ is optically thick to photons having energy greater than 13.6 eV, and 
produces a discontinuity at the hydrogen Lyman limit, i.e. at an observed wavelength of 
912 (1 + z) 'K. These scarcer "Lyman-limit systems" (LLS) are associated with the extended 
gaseous haloes of bright galaxies near the line-of-sight, and have metallicities which appear 
to be similar to that of Lya forest clouds. 

In "damped Lya systems' the H i column is so large (A^hi ^ lO^'^cm"^, comparable to 
the interstellar surface density of spiral galaxies today) that the radiation damping wings of 
the Lya line profile become detectable. While relatively rare, damped systems account for 
most of the neutral hydrogen seen at high redshifts. The typical metallicities are about 10% 
of solar, and do not evolve significantly over a redshift interval 0.5 < z < 4 during which 
most of today's stars were actually formed. 

Except at the highest column densities, discrete absorbers are inferred to be strongly 
photoionized. From quasar absorption studies we also know that neutral hydrogen accounts 
for only a small fraction, ~ 10%, of the nucleosynthetic baryons at early epochs. 



DISTRIBUTION OF COLUMN DENSITIES AND EVOLUTION 

The bivariate distribution /(iVHi,-2) of H i column densities and redshifts is defined by the prob- 
ability dP that a line-of-sight intersects a cloud with column density iVni in the range (iAni, at 
redshift z in the range dz, 

dP = /{Nni, z)dNmdz. (15) 

As a function of column, a single power-law with slope —1.5 appears to provide at high redshift 
a surprisingly good description over 9 decades in Ahi. i e- from 10^^ to 10^^ cm^^. It is a 
reasonable approximation to use for the distribution of absorbers along the line-of-sight: 

f{Nni,z)=AN^tH^ + ^y- (16) 

Lya forest clouds and Lyman-limit systems appear to evolve at slightly different rates, with 
7 = 1.5 ± 0.4 for the LLS and 7 = 2.8 ± 0.7 for the forest lines. Let us assume, for simplicity, a 
single redshift exponent, 7 = 2, for the entire range in column densities. In the power-law model 
(p^) the number N of absorbers with columns greater than Ahi per unit increment of redshift is 

-r = / /(Ki' = 2AN^^-'{1 + z)\ (17) 

az JNui 
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Fig. 1. — High resolution (A/AA = 37,0000) spectrum of the Ze = 3.50 quasar Q1159+123, taken with 
the Keck High Resolution Spectrograph (exposure time 8 h). The data are taken from Songaila (1998). The 
Lyo! forest is clearly seen in absorption blueward of the atomic hydrogen Lya emission line from the quasar 
(the broad peak at 5470 A), and is produced by resonant Lyo; scattering in gas clouds along the line-of-sight 
between us and the quasar. The spectrum shows a Lyman-limit system just shortward of 4150 A, i.e. close 
to the quasar emission redshift. Most of the features between the Lya and C iv emission (the other broad 
peak just below 7000 A) are C iv intergalactic absorption lines. 
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A normalization value A — 4.0 x 10^ produces then ~ 3 LLS per unit redshift at 2; = 3, and, 
at the same epoch, ~ 150 forest lines above A^hi = lO^^-^cm"^, in reasonable agreement with 
the observations. 



If absorbers at a given surface density are conserved, with fixed comoving space number 
density n — nQ{l + z)^ and geometric cross-section E, then the intersection probabihty per 
unit redshift interval is 

^ = E„- = E„„(l + .)3-. (18) 

If the Universe is cosmologically fiat, the expansion rate at early epochs is close to the 
Einstein-de Sitter limit, and the redshift distribution for conserved clouds is predicted to be 

^oc(l + ^f^(x (1 + ^)^/2. (19) 
dz dz 

The rate of increase of /(A^hi, z) with z in both the Lya forest and LLS is considerably faster 
than this, indicating rapid evolution. The mean proper distance between absorbers along 
the line-of-sight with columns greater than A^hi is 



de dz ciV^^ 
dzdN ~ Hon]!i'2A{l + z)^-^' 



For clouds with A^hi > 10^^ cm ^, this amounts to L ~ 0.7 h ^QjJ'''^ Mpc at ^ = 3. At the 
same epoch, the mean proper distance between LLS is L ~ SOh-^n]^^^ Mpc. 



INTERGALACTIC CONTINUUM OPACITY 

Even if the bulk of the baryons in the Universe are fairly well ionized at all redshifts z ^ 5, 
the residual neutral hydrogen still present in the Lya forest clouds and Lyman-limit systems 
significantly attenuates the ionizing flux from cosmological distant sources. To quantify the 
degree of attenuation we have to introduce the concept of an effective continuum optical depth 
Teff along the line-of-sight to redshift z, 

(e-0=e-^««, (21) 

where the average is taken over all lines-of-sight. Negleting absorption due to helium, if we 
characterize the Lya forest clouds and LLS as a random distribution of absorbers in column 
density and redshift space, then the effective continuum optical depth of a clumpy IGM at the 
observed frequency Ug for an observer at redshift Zo is 

Tesii^o, Zo, z) = dz' dNmf{Nm,z){l-e-^). (22) 

Jzo Jo 
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where r = A^hiC"h('^) is the hydrogen Lyman continuum optical depth through an individual 
cloud at frequency v = z/o(l + z)/(l + Zo). This formula can be easily understood if we consider 
a situation in which all absorbers have the same optical depth tq independent of redshift, and 
the mean number of systems along the path is /S.N = J dzdN/dz. In this case the Poissonian 
probability of encountering a total optical depth /ctq along the line-of-sight (with k integer) is 
p(A;ro) = e-^^AA^V(roA;!), and (e^^) = e-^XA;ro) = exp[-AiV(l - e^^«)]- 

If we extrapolate the A^hi^ ^ power-law in equation (|I6|) to very small and large columns, the 
effective optical depth becomes an analytical function of redshift and wavelength, 



reff(i/o, z) = -^lA [^] (1 + zof-'' (1 + zf-'' - (1 + zof-'' • (23) 

Due to the rapid increase with lookback time of the number of absorbers, the mean free path 
of photons at 912 A becomes so small beyond a redshift of 2 that the radiation field is largely 
'local'. Expanding equation ( p3D around z, one gets Tcs{yL) ~ 0.36(1 + 2;)^A2;. This means that 
at z = 3, for example, the mean free path for a photon near threshold is only Az = 0.18, and 
sources of ionizing radiation at higher redshifts are severely attenuated. 
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